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ABSTRACT 

We present H- and ii" s -band imaging of three fields at the centre of 30 Doradus in the 
Large Magellanic Cloud, obtained as part of the Science Demonstration programme 
with the Multi-conjugate Adaptive optics Demonstrator (MAD) at the Very Large 
Telescope. Strehl ratios of 15-30% were achieved in the Tig-band, yielding near-infrared 
images of this dense and complex region at unprecedented angular resolution at these 
wavelengths. 

The MAD data are used to construct a near-infrared luminosity profile for R136, 
the cluster at the core of 30 Dor. Using cluster profiles of the form used by |Elson et al.[ 
we find the surface brightness can be fit by a relatively shallow power-law function 
(7 ~ 1.5-1.7) over the full extent of the MAD data, which extends to a radius of ~40" 
(~10pc). We do not see compelling evidence for a break in the luminosity profile as 
seen in optical data in the literature, arguing that cluster asymmetries are the domi- 
nant source, although extinction effects and stars from nearby triggered star-formation 
likely also contribute. These results highlight the need to consider cluster asymmetries 
and multiple spatial components in interpretation of the luminosity profiles of distant 
unresolved clusters. 

We also investigate seven candidate young stellar objects reported by Gruendl & 
Chu from Spitzer observations, six of which have apparent counterparts in the MAD 
images. The most interesting of these (GC09: 053839.24 —690552.3) appears related to 
a striking bow-shock-like feature, orientated away from both R136 and the Wolf-Rayet 
star Brey 75, at distances of 19'.'5 and 8" (4.7 and 1.9 pc in projection), respectively. 

Key words: instrumentation: adaptive optics - techniques: high angular resolution 
- open clusters and associations: individual: 30 Doradus - Magellanic Clouds 



1 INTRODUCTION 

30 Doradus in the Large Magellanic Cloud (LMC) is the 
largest star-forming region in the Local Group. As an 
archetypal, small-scale 'starburst', 30 Dor provides us with 
an excellent laboratory in which to study star formation and 
stellar evolution, while also giving us potential insights into 
the nature of distant super-star-clusters for which we only 
have integrated properties. Extensive ground-based optical 
imaging and spectroscopy has been used to study the ini- 
tial mass function (IMF), reddening, star- formation history, 



stellar content and kinematics in 30 Dor (e.g. Melnick 1985 
Parker||T993l |Parker fc Garmany||l993l |Walborn fc Blades| 



1997 



Selman 



1998||Bosch et al.|2001p Meanwhile, near-IR 



observations with the Hubble Space Telescope (HST) have 
provided evidence of triggered star-formation, showing the 
region to be a two-stage starburst ( jWalborn et al.||1999| 
[2502| . 

At the centre of 30 Dor is the dense star cluster R136, 
with stellar ages for the most massive stars in the range of 
1-2 Myr ( Massey fc Hunter|1998 1 and a total stellar mass in 
the range of ~0.35-lx 10 D M Q (e.g. [Mackey fc Gilmore|2003 



|Noyola fc Gebhardt|2007| |Andersen et al.||2009| >, depending 
on the low-mass form of the mass function, putting it on a 
par with some of the clusters found in starburst and interact- 
ing galaxies such as M51, M82 and the Antennae. However, 
the core of R136 is too dense for traditional (seeing-limited) 
ground-based techniques. 
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Only with the arrival of HST was R136 resolved in op- 
tical and UV images ( Campbell et al.|1992 de Marchi et al.| 



CAO, Amorim et al. 20061, which operates over the J, 



1993] |Hunter et aT7||1995[ |1997| |Sirianni et al.||2000[ ), with 



follow-up spectroscopy revealing a hitherto unprecedented 
concentration of the earliest O-type stars ( |Massey fe Hunter 
|1998[ ). Star counts from the optical HST images revealed 
that the luminosity profile of R136 appears to be best de- 
scribed by two components, with a break at 10" (Mackey & 
|Gilmore|[2003l ). New results from F160W imaging (equiva- 
lent to -ff-band) with the HST Near Infrared Camera and 
Multi-Object Spectr ometer (NICMOS) fit the profile with 
a single component ( Andersen et al.||2009 l, but only in the 
inner 2pc (8'.'25). Whether the second component seen in 
the optical data is a manifestation of the 'excess light' pre- 
dicted to originate from rapid gas removal in the early stages 
of cluster evolution (Bas tian fe Goodwin||2006 1 remains an 
open question as there is significant and variable extinction 
across the cluster. 

As part of the technology development plan towards the 
European Extremely Large Telescope (E-ELT), the Multi- 
conjugate Adaptive optics Demonstrator (MAD, Marchetti| 



et al. 2007[ ) was developed as a visiting instrument for the 



Very Large Telescope (VLT). The offer of Science Demon- 
stration (SD) observations with MAD presented the perfect 
opportunity to obtain imaging of R136 at unprecedented an- 
gular resolution at near-IR wavelengths. In particular, MAD 
has the power to penetrate the gas and dust more success- 
fully than in the optical HST images, at comparable angu- 
lar resolution, to provide empirical constraints on the outer 
component on the luminosity profile. Determining if R136 
is an expanding group or a dynamically-stable star cluster 
would serve as an important ingredient in the debate on the 



importance of 'infant mortality' of young clusters (e.g. Lada 



fc Lada|2003]|Bastian et al.|2005| [Fall et al.|2005||Goodwin 
fc**FiastTeuipfJ06nBamnga fc Kroupa 2007; Bastian et al 



2008) 



Here we present the MAD SD observations, which de- 
liver a cleaner point spread function (PSF) than NICMOS, 
at finer angular resolution, and over a larger total field. The 
instrumental performance of MAD is discussed in Section [2] 
with the astrometric and photometric calibration of the data 
detailed in Sections [3] and [4] respectively. Following discus- 
sion of seven of the candidate young stellar objects (YSOs) 
reported by |Gruendl fc Chu| ( [2fJ09] hereafter GC09) that lie 
within the MAD fields (Section |5| , we then investigate the 
radial luminosity profile of R136 via a combination of star 
counts and integrated-light measurements (Section [6| , dis- 
cussing its implications in the context of cluster formation 
and observations of distant unresolved clusters. 



2 OBSERVATIONS & DATA REDUCTION 

MAD employs three Shack-Hartmann wavefront sensors to 
observe three natural guide stars (NGS) across a 2' circu- 
lar field, thereby allowing tomography of the atmospheric 
turbulence. The turbulence is then corrected using two de- 
formable mirrors (operating at ~400 Hz), one conjugated to 
the ground-layer (i.e. km), the second conjugated to 8.5 km 
above the telescope. 

The high resolution, near-IR camera used with MAD is 
the CAmera for Multi Conjugate Adaptive Optics (CAM- 



H, and K B bands, with critical (2 pixel) sampling of the 
diffraction-limited PSF at 2.2/xm. The detector is a 2kx2k 
Hawaii-2 HgCdTe array, with a pixel scale of 0'.'028/pixel, 
giving a field-of-view of 57"x57". A useful feature is that 
the camera can be moved within the 2' field without re- 
quiring positional offsets of the telescope, meaning that the 
adaptive optics (AO) loop can remain closed. 

The MAD data presented here were obtained from one 
of twelve SD programmes, observed in November 2007 and 
January 2008. H- and A's-band images were obtained for 
three pointings in the inner region of 30 Dor, as shown in 
Figure [I] The central co-ordinates for the CAMCAO ob- 
servations of Field 1 were a = 05 h 38 m 46.5 s , 5 = -69°05'52" 
(J2000.0). Fields 2 and 3 were offset from this first point- 
ing by approximately —25" in right ascension and ±25" in 
declination. The three NGS used for wavefront sensing are 
also shown in Figure [I] these are Parker #952 (GS1), #499 
(GS2), and #1788 (GS3 ) with V= 12.0, 11.9, and 12.0, re- 
spectively (Parker 1993). The combined 7f a -band image of 
Field 3 is shown in Figure [2] 

The observations are summarised in Table [TJ The de- 
tector integration time (DIT) for all of the observations was 
2 s, with 30 integrations (NDIT) for each exposure. Batches 
of three (H) and six (K s ) object and sky frames were inter- 
leaved in an A-B-A-B-A-B-A pattern, yielding total expo- 
sures of 12min for each field in the H band, and 24min in 
K s . Each science exposure within each batch was dithered 
by 5" - although this reduces the effective area of the final 
combined images, it was intended to minimise the impact of 
bad pixels and cosmetic features from the array. Given the 
vast spatial extent of 30 Dor, the sky offsets were somewhat 
large (+12s of right ascension, +13' in declination) to en- 
sure they were uncontaminated by nebulosity. Observations 
were halted mid-way through the K s exposures for Field 
1 on 2008 January 7 due to bad weather, but were com- 
pleted the following night. Note that the LMC never rises 
above an altitude of approx. 45° as viewed from Paranal, i.e. 
the minimum zenith distance of the observations was ~45°, 
providing a good test of the MAD performances at moder- 
ately large airmass. Indeed, the airmass of the observations 
ranged from 1.4 to 1.6. The range of seeing values for each 
field, as measured by the Differential Image Motion Monitor 
(DIMM) at Paranal, are given in Table [T] 

The MAD data were reduced with standard iraf rou- 
tines, using calibration frames from the SD runs to correct 
for the dark current, to flat-field all of the object and sky 
exposures, and to reject bad pixels and cosmic rays. Me- 
dian sky frames were created for each batch of science ob- 
servations using the sky frames observed immediately before 
and/or after, although (in general) the sky background did 
not vary strongly over each sequence of observations. The 
sky-subtracted frames were then aligned with each other and 
combined. At this stage we omitted one or two images with 
significantly poorer image quality in some fields, hence the 
final exposure times in Table [l] Note that no objects were 
saturated in any of the individual science frames. 

2.1 Image Quality & Performance Analysis 

Moderately bright stars were used to investigate the image 
quality in the three fields. The range and mean of the full- 
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Figure 1. V-band image of the central part of 30 Dor (approx 3f0 by 2.'5) from the Wide-Field Imager (WFI) on the ESO/Max Planck 
Gesellschaft 2.2-m telescope (from observations by J. Alves, B. Vandame & Y. Beletsky). North is towards the top of the figure, east 
towards the left. The reference stars used for the MCAO correction are shown (GS1, 2 & 3), together with the spatial extent of the 
dithered H-band observations of the three fields. 




Figure 2. MAD K s -b&nd image of Field 3. Compared to the overlay on Figure^the dithered regions are cropped, with a final field-of-view 
of ~52" X 46". 
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Table 1. Summary of the VLT-MAD observations in 30 Doradus. The total exposure times quoted are for the final combined images. 



Pointing 


Band 


Date 


Total Exp. 


DIMM range 


Image FWHM 


<FWHM> 








[min] 


["] 


["] 


["] 


Field 1 


K a 


2008/01/07 & 08 


22 


0.4-1.8 


0.10-0.13 


0.11 


Field 2 


K s 


2008/01/07 


24 


0.5-1.1 


0.08-0.10 


0.09 


Field 3 


K s 


2007/11/27 


23 


0.6-1.0 


0.10-0.20 


0.14 


Field 1 


H 


2008/01/08 


12 


0.3-0.6 


0.10-0.12 


0.11 


Field 2 


H 


2008/01/08 


12 


0.9-1.1 


0.08-0.11 


0.09 


Field 3 


H 


2008/01/08 


11 


0.6-1.6 


0.08-0.15 


0.12 



width half maxima (FWHM) obtained from 32 stars evenly 
distributed across each co-added image is summarised in Ta- 
ble [l] Maps of the Strehl ratio in the K s -band images are 
shown in Figure [3] More detailed Strehl maps (showing the 
relative positions of the NGS) and FWHM maps in both 



photometric bands were presented by Campbell et al. ( 200811 . 

The best Strehl ratio achieved in the 7f s -band (25-30%) 
was in the regions closest to the NGS (as one would expect) 
in Fields 2 and 3. The performance in Field 2 is particularly 
good, with an average FWHM of 0'.'09 in both the H and 
K s images (compared to diffraction limits [A/D] of approx 
0'.'04 and 0'.'06, respectively). Although the Strehl is lower in 
Field 1, this was in the best position with respect to all three 
NGS and the correction is very uniform across the combined 
~50"x50" image. This is in strong contrast to Field 3 in 
which the performance is very good in one corner, but then 
steeply declines away from the NGS - more in keeping with 
'classical' AO observations. In general, the performances are 
comparable with those found from other MAD observations 
(e.g. |Marchetti et al.|2007| |Bouy et al.|2008"l >. 



3 ASTROMETRIC CALIBRATION 

Astrometric calibration of each field was undertaken us- 



ing the Selman ( 1998 \ catalogue, recently recalibrated by 
Brian SkifQ These positions are dependent on the preci- 
sion of the 2MASS/UCAC2 positions, combined with the 
plate-scale/seeing of the original New Technology Telescope 
(NTT) observations. The quoted precision on the new Skiff 
astrometry is ~0'.'l (cf. the 0'.'028/pixel delivered by MAD). 
Visual matches of stars between the Skiff catalogue and the 
MAD fields were used to define ~40 well distributed astro- 
metric standards in Fields 1 and 2. The astrometric calibra- 
tion of Field 3, in which there are greater PSF variations, 
was achieved using 60 visually-matched stars. 



4 PHOTOMETRIC CALIBRATION 

Photometry was performed using standard daophot rou- 
tines in iraf on the combined images. Due to the variations 
in the PSF delivered by the AO system across the field, PSF- 
fitting photometry was used, with the model PSF (penny2) 
allowed to vary quadratically across the fields. Roughly 30 
bright, isolated stars were used to create a model PSF for 



each combined image, with neighbouring stars subtracted if 
they were contaminating the selected star. A fitting radius of 
(2 x FWHM)— 1 pixels was found to provide the best fits. Ob- 
jects with 1.5 < sharpness < —1.5, or with x < 4 (where \ is 
the root-mean-square of the residuals that remain) were re- 
jected from the final catalogues to reject non-stellar objects, 
residual cosmic rays etc., as were objects with instrumental 
magnitude errors greater than 0.25 m . 



4.1 Zero-point Calibrations 

4.1.1 MAD vs. 2MASS 

Sources from the Two Micron All Sky Survey (2MASS, 
|Skrutskie etwST1|2l)06"| ) were identified in the MAD frames 
for photometric calibration. However, at the excellent angu- 
lar resolution delivered by MAD many of the 2MASS 'stars' 
are resolved into asymmetric sources or multiple components 
(e.g. Fig. 3, Momany et al. 2008 1 . As such, only isolated, 



apparently single stars with 2MASS photometric qualities of 
either 'A', 'B', or 'C in the relevant band were considered 
for calibration. This limited the number of 2MASS stars in 
each MAD field to eight or fewer for calibration. 

Photometric zero points (ZPs) were determined (prin- 
cipally using 'A'-rated 2MASS stars), with stars at the 
bright/faint ends of the MAD images omitted to exclude 
strongly-exposed stars in the MAD frames and poorly- 
exposed stars in 2MASS. The resulting ZPs are summarised 
in Table [2] we are somewhat at the mercy of small number 
statistics, not to mention the issues of stellar density and 
complex nebular emission in and around R136. 

4.1.2 MAD vs. HAWK-I 

To investigate the quality of the calibrations we employed 
i\s-band commissioning data from the VLT High Acuity 



Wide-field K-band Imager (HAWK-I, Pirard et al. 2004 



Casali et al. 20061. HAWK-I is a wide-field, near-IR cam- 
era with four 2k x 2k Hawaii-2 detectors, covering a 7!5 
square field at a pixel scale of 0'.'106. Commissioning images 
of 30 Dor were obtained on 2007 October 3, using Y , J, 
K a and Br7 filters. The HAWK-I exposures were dithered 
around the centre of 30 Dor to build-up a complete mosaic 
for the regiorj^] 

One A" s -band HAWK-I array with a reasonable spatial 
overlap was selected for each MAD image. These frames were 



1 http://cdsarc.u-strasbg.fr /viz-bin/Cat?J/A+A/341/9 



See p. 22 of the December 2007 edition of The Messenger. 
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Figure 3. i\g-band Strehl (%) maps for Fields 1, 2 and 3 (left to right, respectively). North is at the top, with east at the left. 



Table 2. Photometric zero-points obtained from the combined image method. Results are given for the iC s -band images using 2MASS 
sources within the MAD frames, and using the HAWK-I commissioning data, in which the result for Field 2 was bootstrapped using 
overlapping stars with Field 1. The H-band results were obtained using 2MASS sources for Fields 1 and 3, with Field 2 again calibrated 
using overlapping stars with Field 1. 



Band 


Calibration 


Field 1 


Field 2 


Field 3 


K B 


2MASS 


26.78 ± 0.13 


26.61 ± 0.45 


26.95 ± 0.25 


K s 


HAWK-I 


26.69 ± 0.08 


27.04 ± 0.09 


27.07 ± 0.11 


H 


2MASS 


27.09 ± 0.14 


27.27 ± 0.15 


26.93 ± 0.08 



then reduced and analysed using the same methods as for 
the MAD images, i.e. using PSF-fitting photometry. There 
were ~50 stars from 2MASS with 'AAA' JHK a quality rat- 
ings in each of the HAWK-I frames. The resulting K a ZPs 
were applied to all of the objects found in the image and a 
catalogue was created. A comparison of the HAWK-I mag- 
nitudes with those from 2MASS for the calibration stars is 
shown in Figure |4] with dispersions of ±0.12, 0.16, and 0.21 m 
for Fields 1, 2 and 3, respectively. 

These results provided a well-defined ZP between 
HAWK-I and 2MASS, which was then used to bootstrap 
the MAD ZPs from stars overlapping in the images. By com- 
paring the instrumental magnitudes of 30 stars in the MAD 
frames with their corresponding HAWK-I values, new ZPs 
were calculated for each K s MAD field. Note that the stars 
used were hand-picked to be isolated and well within the dy- 
namic range of the MAD observations - using fainter sources 
from the MAD data would have pushed toward the sensi- 
tivity limits of the HAWK-I frames rather than providing 
improved calibration of the MAD data. The internal agree- 
ment of overlapping stars within the three HAWK-I frames 
used for calibration of the three MAD fields was found to be 
better than 0.1 m . 

Stars that appear in the overlap regions of the MAD 
fields were used as a double-check of the ZPs from the 
HAWK-I data. Fields 1 and 3 were in excellent agreement, 
while the ZP for Field 2 gave magnitudes that were 0.07 m 
brighter, although still within the calculated uncertainties; 
the ZP for Field 2 was corrected accordingly, resulting in 
i\s-band magnitudes on a common system. 

Unfortunately there were no matching HAWK-I ob- 
servations in the ff-band to employ a similar calibration 




s HAWK-I 



Figure 4. Residuals of the calibrated HAWK-I Kg-band photom- 
etry compared with the 2MASS values. Field 1: open squares; 
Field 2: open triangles; Field 3: crosses. 



method. Informed by the behaviour of the if s -band results, 
we adopt the 2MASS ZPs from Table [2] for Fields 1 and 
3, which have much smaller standard deviations than found 
for Field 2. Instrumental magnitudes were then compared 
for stars overlapping between Fields 1 and 3 as done for K a ; 
the resulting offsets were in excellent agreement with those 
expected from the ZPs from Table [5] (i.e., AH~0.16 m ). We 
then compared instrumental magnitudes for stars overlap- 
ping between Fields 1 and 2. These offsets were then used 
to calculate a new mean ZP for Field 2, which is within the 
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Table 3. Photometric zero-points obtained from the individual 
frame method. As in Table [2] the A" s -band results are derived 
from HAWK-I, with the //-band results from 2MASS. In both 
cases Field 2 was calibrated using stars overlapping with Field 1. 



Band 


Field 1 


Field 2 


Field 3 


K B 


23.40 ± 0.07 


23.62 ± 0.08 


23.59 ± 0.11 


H 


24.45 ± 0.14 


24.37 ± 0.14 


24.33 ± 0.08 



quoted uncertainties from the 2MASS comparison, but has 
a much reduced error. 



4-1-3 PSF Fitting on Individual Frames 

If the PSF of a star varies subtly between different frames (as 
is likely with AO-corrected imaging) the resulting combined 
PSF may broaden slightly and be noisier. To investigate if 
this effect was the dominant source of our ZP uncertain- 
ties, we employed a second photometric method to see if we 
could improve the calibrations. We used PSF-fitting subtrac- 
tions on the individual frames, then used Dr. Peter Stetson's 
DAOMATCH and daomaster packages to create a mean cat- 
alogue. This method potentially provides more precise PSF 
fitting, leading to improved fidelity of the final photometric 
catalogue. It also allows the user to set the number of frames 
in which a source must be detected for its inclusion in the 
final catalogue. This has the big advantage that objects in 
the dithered regions (previously discarded when the com- 
bined image is subset) can now be included, increasing the 
spatial extent of the final source catalogue, and also leads to 
fewer rejections of stars affected by cosmics or bad pixels in 
only one frame. 

The uncertainties on the ZPs obtained from analysis of 
the individual frames were (effectively) unchanged compared 
to those found from our previous best efforts (see Table [3] cf. 
Table [5| . This suggests that the uncertainties in the mag- 
nitudes of the calibration stars are dominanting those from 
the photmetric methods employed. A more tangible gain is 
the increased number of sources detected due to the reten- 
tion of the dithered regions not observed in all frames (albeit 
with lower sensitivity) . Following application of the ZP cali- 
brations, a comparison between the combined- frame and in- 
dividual catalogues was undertaken for each pointing/band, 
finding mean offsets less than 0.05 m in all six instances. 

The ZPs adopted for calibration of the MAD fields for 
the final catalogue, created from the individual frames, are 
listed in Table [3] In summary, the _K" a -band ZPs were cali- 
brated using HAWK-I, with Field 2 bootstrapped using over- 
lapping stars in Field 1. The //-band ZPs were calibrated 
using 2MASS, also with Field 2 adjusted from comparisons 
with Field 1. Figure [5] compares the final MAD magnitudes 
with those of the HAWK-I calibration stars, with the dis- 
persions around the mean quoted in Table [3] 

As an external check on the MAD magnitudes we com- 
pared our results for the 12 (visually-matched) overlapping 
stars in the NICMOS observations from IBrandner et al.l 
( |2001| . Nine of these stars were in Field 2 (with two only 
observed in the J^-band due to the dither pattern), with 
the remaning three in Field 3. Accounting for the slight 
differences in the photometric filters (the HST data were 




MAD 



Figure 5. Residuals of the calibrated MAD Kg-band photometry 
from the individual frame method compared with the HAWK-I 
values for hand-picked reference stars (see text for details) 



transformed to the CTIO photometric system), the mean 
differences (MAD — HST) and standard errors are in excel- 
lent agreement AH = -0.05 ± 0.03 and AJs" s = -0.04 ± 0.05 
(std. err.). 

4.2 Colour-magnitude diagrams 

Colour-magnitude diagrams (CMDs) were created from cat- 
alogues from both photometric methods (i.e. from the com- 
bined and the individual frames), excluding objects within a 
radius of 2'.'8 of the centre of R136, where crowding/blending 
become dominant. 

The individual frame method yields an extra 600-1,000 
objects in the CMDs compared to the analysis of the com- 
bined frames, mainly due to the increase in the effective 
field-of-view. When CMDs from the individual method are 
compared over the same spatial region as the combined 
CMD, Fields 2 and 3 contain slightly fewer objects, pre- 
dominantly at fainter magnitudes, suggesting the individual 
frame method results in slightly diminished sensitivity. In 
contrast, Field 1 still contains fewer objects in the combined 
CMD, which is a legacy of a light leak down one side of 
the images that affected some of the MAD SD observations. 
Field 1 was the most affected by this light leak in our ob- 
servations, with the image quality parameters in daophot 
rejecting more sources in analysis of the combined frame in 
this region than in the individual frames. 

In addition to the inclusion of the dithered regions, the 
individual frame method reduces the random star-to-star 
uncertainties in the final catalogues. CMDs from the individ- 
ual frame method are shown for all three fields in Figure [6] 
The location of the main-sequence is in good agreement for 
all three fields. 

There is significant and variable extinction toward 
R136, with | Andersen et al.| ( |2009[ > adopting a median ex- 
tinction of Av = 1.85 for massive stars (M = 7-20 Mq), cor- 
responding to Ak ~ 0.21 (assuming, for the purposes of an 
order of magnitude estimate of the IR extinction, the stan- 
dard galactic extinction law from Ricke & Lebofsky 1985 1 
and Ah ~ 0.32 (using the scaling from |Indebetouw et al.| 
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20051, yielding a relatively small near-IR extinction term of 
E(H-K S ) = 0.11. Thus, the H-K s colours from MAD will, 
in general, not be differentially reddened by much more than 
the photometric errors. CMDs were created for different ra- 
dial bins using the combined catalogue from all three fields 
and these show no significant offset in the locus of the main 
sequence. 

For Field 1 in Figure || we include the unreddened 
logi = 3.0 isochrone (effectively illustrative of the zero-age 
main sequence) from Lejeune & Schaerer (20011, adopting 



the tracks with the metallicity relevant to that of the LMC 
(i.e. those from [Schaerer et al.||1993] ). Note that the offset 
from the MAD data is in agreement with the reddening es- 
timated above. 
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4.3 Photometric Completeness 

Completeness tests were undertaken in both bands for all 
three fields. Using the STARLIST and ADDSTAR routines in 
iraf, 100 artificial stars with magnitudes varying from 14 
to 24 were distributed uniformly across each combined (and 
subset) image. As before, PSF-fitting was performed, using 
the same settings and a PENNY2 model PSF which is al- 
lowed to vary across the field. Using the known positions 
of the added stars, the number recovered from the images 
within a two-pixel search radius was found. If more than one 
object was found within that small radius, the object with 
the smallest magnitude difference was considered the arti- 
ficial star. Stars with magnitude differences of greater than 
±0.5 m were cut from the detected list, with the intention of 
excluding mistaken matches, or artificial stars which have 
been superimposed on genuine stars within the image (thus 
increasing its magnitude). 

These tests were repeated 1,000 times, until 100,000 ar- 
tificial stars had been added to each image. This gives a 
ratio of ~30 to 50 artificial stars for each observed object. 
The ratio of artificial stars introduced to those detected was 
measured for each image, as shown in Figure [7] Field 2, 
which does not include R136 in the combined image, has a 
slightly more uniform completeness at brighter magnitudes 
while going slightly deeper in both bands. 

The 50% completeness level in the _ftT s -band frames is 
K s = 19.45 in Field 1, 19.85 in Field 2, and 18.90 in Field 3, 
corresponding to initial main-sequence masses of approxi- 
mately 3.5, 2.8, and 4.7 M , respectively, when compared 
to the youngest isochrone from |Lejeune fc Schaerer|p00l] >. 

For Fields 1 and 3 the completeness was also found for 
both bands as a function of radius from the cluster (adopt- 
ing 100 pixel radial bins), an example of which is shown in 
Figure [8] These radial completeness profiles are taken into 
account in construction of the radial luminosity profile in 
Section [6] 



5 SPITZER YSO CANDIDATES 

From observations with the Spitzer Space Telescope, GC09 
presented a catalogue of potential young stellar objects 
(YSOs) in the LMC. Prior to investigation of the luminos- 
ity profile of R136, we first investigated the seven candidate 
YSOs from their catalogue that are within the MAD fields. 
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Figure 6. CMDs from the mean catalogues created from the in- 
dividual frames method, excluding the central 2'.'8 of R136. The 
dashed horizontal lines mark the 50% completeness level (cf. Fig- 
ure [TJl in the common regions included in all positional dithers. 
Overplotted in the Field 1 panel is the unreddened youngest 
isochrone (logt = 3.0) from Lejeune & Schaerer (2001), with the 
marked points indicating, from bright to faint magnitudes, initial 
masses of 60, 50, 25, 15, 9, 5, 3, and 2M . 



The Spitzer data were primarily obtained under the aus- 
pices of the Surveying the Agents of a Galaxy's Evolution 
(SAGE) LMC survey ( jMeixner et al.| |2006), but were also 
supplemented by smaller, more targeted programmes from 
the Spitzer archive. 

Following initial colour cuts to exclude the majority of 
stars and background galaxies, GC09 investigated the nature 
of the candidate YSOs via inspection of their morphologies 
and spectral energy distributions, including comparisons at 
other wavelengths. This resulted in five categories for the 
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Figure 7. H (upper panel) and K B (lower) completeness profiles 
for Fields 1, 2 and 3, as indicated. 



colour-selected Spitzer sources: (1) evolved stars, (2) plan- 
etary nebulae, (3) background galaxies, (4) diffuse sources, 
and (5) 'definite', 'probable' and 'possible' YSO candidates. 
This is a somewhat different approach to the selection of 
candidate Y SOs from the SAGE LMC surve y by |Whitney 



et al. (2008 see further discussion by GC09). Vaidya et al 
( ]2009| ) have since compared the Spitzer data with archival 
HST Ha images and new (seeing-limited) near-IR imaging, 
enabling more detailed description of the local environments 
around 82 of the YSOs from GC09. 

The positional uncertainty on the Spitzer astrometry 
is approximately ±0'.'2 (Dr R. Gruendl, priv. comm.), but 
in the crowded environment of 30 Dor, one might expect 
that blended sources and nearby gas might lead to greater 
uncertainties. Indeed, the typical angular resolution of the 
Spitzer observations ranged from 1'.'7 to 2'.'0 over the four 
InfraRed Array Camera (IRAC) bands for the SAGE obser- 
vations ( Meixner et al.|2006 1. Thus, it is not surprising that, 
in the case of two of the 'diffuse' sources from GC09, we find 
apparent counterparts in the MAD images offset by a radius 
of approximately 1". 

In light of these positional uncertainties, there remains 
one candidate from GC09 which does not have an obvious 
close counterpart within a 2'.'0 radius: the 'diffuse' candidate 
053842.69 -690623.7. There is a faint source (K s = 18.1, 
H—K s = 0.39) at a distance of ~f/.'75, which is otherwise 



unremarkable; a brighter counterpart can be ruled out. We 
note that of the four Spitzer-IRAC bands, GC09 only re- 
ported magnitudes for this source at 4.5 and 8.0 /im, sug- 
gesting a less robust detection than for the other six sources 
discussed here which were detected in all four IRAC bands. 

The remaining six YSO candidates are now discussed 
in turn, with their locations overlaid on the V-band WFI 
image in Figure [9] to illustrate their positions relative to 
R136. All except 053836.48 - 690524.1 feature in the com- 



parison by Vaidya et al. (20091, with four also within the 
(seeing-limited) near-IR imaging from |Rubio et al.| ( |1998| 
hereafter RBW98). Spitzer spectroscopy of three of these 
sources [053839.24 - 690552.3, 053839. 69 - 690538.1 and 
053845.15 -690507.9) was presented by |Seale et aT^p009) , 
with each classified in their 'PE Group', which are seen 
to display polycyclic aromatic hydrocarbon (PAH) emission 
features combined with fine-structure lines from ionised gas. 



5.1 053836.48 - 690524.1: 

Reported as a 'diffuse' source by GC09, this object is on 
the western edge of the Field 2 images. Subset (5" x 5") 
H and Jf a -band images are shown in Figure |10| in which 
the cross in the i/-band image indicates the position from 
GC09 and the intensity scaling is the same in both bands. 
The circle (centred on the Spitzer position in the if-band 
image) indicates the typical angular resolution of IRAC in 



the SAGE survey, i.e. l'.'7 at 3.6 /xm (Meixner et al. 20061 



Just below the Spitzer position in the i/-band image is a 
small cluster of bad pixels - the likely YSO counterpart is 
to the east and is more clearly seen in the _/f s -band image. 

Apart from the obvious asymmetry in the A^-band im- 
age, we were somewhat cautious to match the MAD object 
to the GC09 source. However, |Kim et al.| ( |2007| also report a 
candidate YSO (their '30Dor-15', also from IRAC observa- 
tions) that is only ~0'/33 from the MAD counterpart, which 
has K s = 16.24, H — K s — 2.080 



5 . 2 053838. 35 -690630.4 ■' 

Classified as a 'diffuse source' by GC09, there is a poten- 
tial counterpart (K s = 14.37 and H-K s = 0.35) in the MAD 



images at a distance of 1" (Figure 111. The source appears 
to have an associated 'plume' of nebulosity which spirals to 
the south-west for at least 0'.'5. This nebulosity can also be 



seen in the NICMOS image from Brandner et al. ( 2001 their 



30Dor-NIC01 frame), but did not feature in their discussion. 

The 2MASS J-band magnitude quoted by GC09 for this 
source is 13.93 (2MASS source: 053838.5-6906297). Over- 
plotting the 2MASS catalogue on the MAD image reveals 
reasonable astrometric agreement (0'.'2) considering the lim- 
ited angular resolution of 2MASS, with upper limits of 
H > 12.88 and K s > 12.05 (i.e. 'U' rated). 

This source also corresponds to IRSW-98 from RBW98 
(which, in turn, is component 'a' of the knot comprising 



3 We note the candidate YSO '30Dor-22' from |Kim et al.|||2007^ 
is, on the basis of their published astrometry, Mk34 l |Melnick| 
[1985}, a well-studied WN star. We are puzzled that Kim et al. 
mention this as one of the sources studied by |Brandnc r et al7| 
|2001[ |, as it does not feature in their NICMOS fields. 
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Figure 8. Kg-band radial completeness for Field 1. Each annulus has a width of 100 pixels (2''8). 'Radius 1' refers to the radius closest 
to the cluster core. 



source W4 from |Rubio et aL][l992|), w ho found K a = 14.02 
and H-K B = 0.45. |Vaidya et al.| ( |200"9| ) note that this YSO is 



in a marginally-resolved H II region, suggesting that the line 
contribution (at least in the if s -band) is likely contributing 
to the brighter magnitude from RBW98 compared to the 
MAD results. 



5.3 053839.24 - 690552.3: 

The most likely YSO counterpart is at a = 05 h 38 m 39*28, 
5 = -69° 05' 52'.'60 (J2000, ~0'.'45 from the Spitzer posi- 
tion) and is the most immedi ately impressive counter- 
part in the MAD images (Figure [l2 ^1 We find K s = 14.71 
and H—K s = 0.42 for the star, as compared to K B — 13.91, 
H-K s = 0.61 from RBW98 (their source IRSW-118). The 



images from RBW98 were taken in typical seeing of 1" so 
the two sources visible a the centre of Figure [12] will be 
strongly blended, leading to the brighter magnitude. 

Particularly striking is the apparent bow-shock, most 
easily seen in the /f s -band image. This is almost, but not 

4 The source is on the edge of the region common to all H- and 
_K" s -band images, hence the apparent 'banding' where the sensi- 
tivity differs due to the science-frame dithers. 



quite, aligned with the core of R136 approximately 19'.'5 
away. Also of note in the same direction is R134/Brey 75 



I Breysacher 



Smith 



19811, classified as WN6(h) by Crowther fc 



(199711, at a distance of only 8". To date, nothing 
is known regarding the spectral type of the two bright stars 
at the centre of the image (separated by ~0'.'33). This region 
clearly warrants more detailed investigation in the context 
of small-scale triggered star formation. 



5.4 053839.69 -690538.1: 

Classified as a 'definite YSO' by GC09, this object appears 
as a somewhat extended or embedded source in the //-band 
MAD image (Figure |l3), with K s = 13.75 and H-K s = 3.37. 
This source corresponds to IRSW-127 from RBW98 (W9 



from Rubio et al. 1992), for which they found K s = 13.91 



and H-K s = 2.85. 

The bright star l'.'2 to the south is P733/S206, for which 
K s — 14.40 and H—K s = 0.16, suggesting it as a massive star 

Plez||2006] |. The 



(cf. the intrinsic colours from [Martins 



strong contrast in the colours of the two stars is immedi- 
ately obvious from the composite-colour image from RBW98 
(their Figure 1), although they appear slightly blended. 
Blending could perhaps account for the photometric differ- 
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Figure 9. Location of the six YSOs from |Gruendl fc Chu| ( |2009} 
with counterparts in the MAD images, overlaid on the V-band 
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Figure 10. 053836.48 -690524.1: H- and if s -band 5"x5" images 
(left- and right-hand panels, respectively), with north at the top, 
east at the left. The black cross marks the location of the Spitzer 
position for the YSO candidate, with the black circle representing 
the angular resolution of the IRAC data (1'.'7 at 3.6 fim). 



ences (which are in the expected direction), but RBW98 
employed PSF-fitting methods so should be relatively ro- 
bust to such effects at this separation - perhaps our efforts 
are somewhat limited by the apparent asymmetric flux to 
the south of the star. 



5 . 5 05384 3. 52 - 690629. 0: 

Classified as a 'possible YSO' by GC09, the nearest spatial 
match is a source with K s — 17.15 and H—K s = 1.01, as indi- 



east) is P1064/S696, with K B = 16.87 and H-K s = 0.18, and 
for which the spectral type is unknown. This YSO source is 
not within the region observed by RBW98 but, interestingly, 



is noted by Vaidya et al. ( 2009 1 as being in a dark cloud and 



as comprising multiple YSOs in the Spitzer PSF, one with an 
optical counterpart in the Hcv HST images, and one without. 

5.6 053845.15 - 690507.9: 

Classified as a 'definite YSO' by GC09, their published po- 
sition is approximately 0'/4 north of P1222/S116/IRSN-101 
( |Parker||l9"9"3l |Selman||1998[ RBW98), which was classified 
as 03-6 V by Walborn & Blades (19971 from ground-based 



spectroscopy, later revised to 09 V(n)p by Wa lborn et al.| 
(20021 from HST observations; 5"x5" images of the region 



are shown in Figure [T5| Note that PI 222 was not used in 
the astrometric calibration of the MAD frames and yet its 
recovered position agrees with Skiff's astrometry to better 
than two pixels. 

This source lies within the dense nebular region referred 
to as 'Knot T by |WalbornHT986] ), highlighted as young mas- 
sive stars just emerging from their natal cocoons and imaged 
at optical and near-IR wavelengths with HST by |Walborn] 
|et al.| ( |1999"l|2002[ ). The relatively bright Spitzer magnitudes 
quoted by GC09 are in keeping with P1222 as the most 
plausible counterpart. Indeed, the spectral energy distribu- 
tion shown in Figure 13 of GC09 reveals increased mag- 
nitudes bluewards of 1 /im, as well as redwards, consistent 
with a massive star with an IR excess. From the MAD im- 
ages we find K s — 14.31 and H—K s = 0.13, as compared to 
K s = 14.45 and H-K s = 0.03 from the ffST-N ICMOS ob- 
servations of Knot 1 by Brandner et al.| ( |2001[ ) . 

|Vaidya et al.| ( |2009 1 note the source as residing in a 
resolved H II region and, on the basis of the Ha luminosity 
(uncorrected for extinction) provide a lower bound on the 
inferred spectral type of the YSO of B0. 5; evidently from 
the classification of Wa lborn fc Blades) the star is slightly 
hotter. 



5.7 Spatial distribution and masses of YSOs 

It is not surprising that five of the six YSO candidates dis- 
cussed here lie to the north and west of R136 (Figure [9| 
This is the location of the 'second generation' of massive 
stars ( jWalborn et al.|1999||2002[), with significan t molecular 



cated in Figure 14 The adjacent object (0'.'25 to the north- 



material still remaining (e.g. Werner et al.|1978 Johansson 
|et al.|1998| >. The sixth source (053843.52 - 690629.0) is spa- 
tially distinct from the regions of triggered star formation 
and is the most ambiguous of the YSOs discussed here; with 
Mk ~ —1-4, it is also likely to be the lowest mass object. 

The absolute A" s -band magnitudes of our four bright- 
est YSO counterparts are all bright, with Mk ~ —3.75-4.75, 
suggesting that they are bona fide candidates for massive O- 
or early B-type stars, depending on the degree of obscura- 
tion/reddening. The most intriguing in this context is the 
very red source 053839.69 —6905381 which, given the un- 
certainties of the Spitzer astrometry, would certainly benefit 
from ground-based spectroscopy to elucidate its nature and 
its relationship with P733 nearby. |Rubio et al.| ( |2009[ ) have 
recently reported a dense molecular cloud toward this source 
(their IRSW-127) suggesting that it may well be a massive 
star still heavily enshrouded in its natal cocoon. 
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6 RADIAL LUMINOSITY PROFILE OF R136 

We now compare the radial luminosity profile for R136 
from the MAD data, constructed using a combination of 
integrated-light measurements and star counts, with the op- 
tical HST results. 



■ ■* 



Figure 11. 053838.35 -690630.4: 5" X 5" H- and AT s -band im- 
ages. 




Figure 12. 053839.24 -690552.3: H- and if s -band. Note the 
particularly striking bow-shock feature in the K s image. 
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Figure 13. 053839.69 -690538.1: H- and if a -band. 
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Figure 14. 053843.52 -690629.0: H- and i^-band. 




Figure 15. 053845.15 -690507.9: H- and if a -band. 



6.1 Integrated Light Profiles 

Integrated-light measurements were used to investigate the 
innermost part of the radial luminosity profile in Fields 1 
and 3. These sum the flux within annuli around the core, 
where the radius of each annulus is defined such that the 
total flux in each is approximately constant. An idl pro- 
gram was written specifically for this purpose, with standard 
daophot routines in iraf also used as a comparison. The 
idl program calculates the integrated light by re-sampling 
each image as a function of azimuth (<j>) and radius from 
the cluster centre, to enable a smoother, more robust calcu- 
lation at smaller radii than possible with the (albeit small) 
sampling of the pixels - in practice the results are essentially 
the same as those obtained with iraf. These measurements 
adopt a constant sky value, which is defined well away from 
the cluster core. 

To estimate the uncertainties in the integrated-light 
profile, each annulus is split into azimuthal sections and 
the density is calculated in each; the standard deviation of 
the results around each annulus then provide an estimate 
of the statistical uncertainties, provided the cluster is sym- 
metric. However, visual inspection of the images shows that 
the cluster is not symmetric, particularly at larger radii (see 
Section |6.2.1[ ) , suggesting that the variations in the uncer- 
tainties are shaped by both statistical errors and asymmetric 
variations of the cluster. 

The surface brightness profiles were then fit with |Elson"l 
|et al.| |I987| EFF) cluster profiles of the form: 



M(r)= M0 (l+^r 7/2 , 



(1) 



where a is the scale radius that defines where the flat inner 
part (the core) turns into a power-law profile with an index 
of —7. The results are given in Table |4j both in terms of 
angular extent on the sky and in parsecs (adopting a distance 
modulus to the LMC of 18.5). The profile fits are shown in 



Figure 16 in which the lower three profiles are offset (in 
steps of two magnitudes) for clarity. 

There is hardly any signature of a core, with the surface 
brightness profile fit by a power-law function over the full 
range of the data. The scale radius (a) obtained is effectively 
~0.1" (0.025 pc) which is, in practice, the resolution limit 
of these data; i.e. 0.1" provides an upper limit to the core 
radiu^] Both the results for the core and the slope of the 
power-law fit (i.e. 7~ 1.6) are in good agreement with past 
optical HST results ( jCampbell et al.||1992| |Hunter et al.| 
1995} and the new NICMOS results ( jAndersen et al.|2009[ ). 



5 Note that the core radius, r C orc, is usually defined as the radius 
where the surface brightness profile drops to half its central value. 
For an 



Elson et al. 



profile it follows that r c , 
the values of 7 we find here, r core fti 1.5 a. 



:a ( 2 2/7_i). For 
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Table 4. Structural parameters for R136 from EFF fits to the MAD integrated-light profiles. The italicised results are from fits to the 
combined luminosity profiles (see Section |6.3|, For convenience, results are quoted in direct observables and their physical size in pc. 



Field 


/^o [mag arcsec 2 ] 


a[»] 


7 


rc ["] 


H Fl 


6.37±0.12 


0.13±0.03 


1.67±0.15 


0.15 ±0.04 




6.36± 0.12 


0.12± 0.02 


1.57± 0.05 


0.14 ±0.03 


H F3 


6.59±0.11 


0.12 ±0.03 


1.58±0.13 


0.14 ±0.03 




6.59± 0.11 


0.12± 0.02 


1.57± 0.05 


0.13± 0.02 


K Fl 


6.35±0.14 


0.13±0.04 


1.58±0.15 


0.15 ±0.04 




6.35± 0.14 


0.12± 0.02 


1.48± 0.06 


0.15± 0.03 


K F3 


6.39 ±0.10 


0.14±0.03 


1.63±0.14 


0.16 ±0.04 




6.38± 0.11 


0.12± 0.02 


1.56± 0.06 


0.15± 0.03 




Field 


MO [mag pc 2 ] 


a [pc] 


7 


r c [pc] 


H Fl 


3.293±0.121 


0.031 ±0.008 


1.67±0.15 


0.035 ±0.009 




3.282± 0.121 


0.028 ± 0.005 


1.58± 0.05 


0.033 ± 0.006 


H F3 


3.517±0.112 


0.028 ±0.007 


1.58±0.13 


0.034 ±0.008 




3.509± 0.112 


0.021 ± 0.005 


1.57± 0.06 


0.032± 0.006 


K Fl 


3.278 ±0.136 


0.031 ±0.009 


1.58±0.15 


0.037 ±0.010 




3.269 ± 0.136 


0.029 ± 0.006 


1.48± 0.06 


0.035 ± 0.007 


K F3 


3.316 ±0.105 


0.033 ±0.008 


1.63±0.14 


0.039 ±0.010 




3.305 ± 0.105 


0.030 ± 0.006 


1.56± 0.06 


0.036 ± 0.007 



were applied to correct the densities if part of a given anrmlus 
extended past the edge of an image. 

The radial completeness tests (e.g. Figure[8| reveal, not 
unsurprisingly, that the first couple of annuli beyond r = 2'.'8 
are also strongly affected by crowding, and are thus subject 
to greater uncertainties than in the outer regions. The 50% 
completeness level in the third radial bin was adopted as 
the faint limit for construction of the profiles from the star 
counts, corresponding to a faint magnitude limit in each 
field in the range of 18. 5-18. 8 m , equating to a main-sequence 
mass of ~5 Mq . The innermost bins were not included in 
construction of the final combined profile. 

6.2.1 Azvmuthal Density Variations 

To investigate the apparent asymmetry in R136 more quan- 
titatively, we calculated the azimuthal density variation be- 
yond the central 2'.' 8. An idl program was used to split the 
images into angular slices, each containing approximately 
equal numbers of stars (after completeness corrections were 
applied). Similar techniques have been used recently to in- 
vestigate the structure and asymmetries seen in lower-mass 
Galactic clusters (e.g. |Gutermuth et al.|2005[|2009| ). 

Due north from the cluster core was set as 0=0°, in- 
creasing anti-clockwise (increasing towards the east), with 
the outer radius set by the distance from the cluster core to 
the northern edge of Field 1. This enables a near complete 
azimuthal sweep in the combined J/-band images of Fields 1 
and 3 (i.e. not the full extent of those shown in Figure [TJ. 
An iterative routine was used to define an azimuthal slice 
in each image (increasing by 1° each loop), which summed 
the star counts in that slice and corrected for the incom- 
pleteness (down to 50%). Once 200 stars are obtained, the 
routine determines the area covered by each slice, calculates 
the density, and then moves on to define the next slice. 

The azimuthal variation around R136 (with 
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Figure 16. EFF fits to the MAD integrated-light profiles, with 
the K s - and H-band (Field 3) profiles shifted to fainter magni- 
tudes for clarity. 



6.2 Star Counts 

Star counts in annuli around the core were used to calculate 
the radial luminosity profile in the regions beyond r = 2'.' 8. 
Stars were binned such that a similar number of objects 
were included per annulus, as suggested by Mai'z Apcllaniz & 



Ubeda 



I 2005). The stars were binned as a function of radius 



and magnitude and completeness terms were calculated. As 
the core is off-centre in each frame, area correction terms 
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r ma jt = 28'.'0, =1,000 pixels) is clearly seen in Figure 17 
demonstrating that this region is far from symmetric, with 
a local minimum around = 90°. Note that the _ff-band 
densities are larger as a consequence of the deeper images 
compared to the _K" a -band observations (see Figure [7]). The 
exact densities are sensitive to the completeness threshold 
employed, so similar calculations were undertaken for 
completeness cuts of 30, 40, and 60% (with the results for 
40% shown in Figure |18) . While there are variations in 
the exact densities as a function of (j) (P ar tly because the 
different depths yield slices which vary in size in terms of 
<f>), the overall trend remains the same. 

Adopting a distance modulus of 18.5 to the LMC, 
r = 28'.'0 corresponds to a projected radius of 6.8 pc. Similar 
trends in the azimuthal profiles are also seen adopting outer 
radii of 600 and 800 pixels (4.1 and 5.4 pc, respectively). At 
smaller radii incompleteness effects become more significant 
(cf. Figure [8| and it is harder to characterise the asymmetry 
meaningfully. Superficial inspection of Figure 2 from |Massey| 
|fc Hunter] ( |1998| ) suggests similar asymmetries in the inner 
region, with luminous stars appearing to be more numerous 
in the eastern half of the central 4'.'55 x 4" 55 HST image. 
However, we also note that the most luminous, massive stars 
will likely have a different relaxation time to the lower-mass 
population. 



6.3 Combined Radial Profile 

The overlap region between the profiles from the two meth- 
ods was used to convert (and normalise) the star counts into 
surface brightness - a reasonable approach provided there is 
no evidence of mass segregation in the cluster (as claimed 
by |Hunter et al.|[l995[ although Brandl et al. 1996 argued 
to the contrary). The average offset in the overlap region 
was found from interpolation between the two profiles, re- 
sampled to an evenly sampled grid of radial bins to avoid 
a bias toward smaller radii (at which the bins are smaller) 
and excluding the innermost bin from the star counts profile 
(to be less sensitive to completeness corrections). The offset 
was then used to convert the star counts into a magnitude 
density, yielding a combined profile. 

EFF fits to the combined profile are shown in Figures (19] 
and |20[ with the fit parameters given in italics in Table [4] 
Meaningful errors on the star counts are hard to estimate 
due to the difficulty of attributing uncertatinties to the com- 
pleteness corrections (i.e. blending/crowding effects) and the 
area corrections (i.e. cluster asymmetries), so we adopt a 
conservative error of ±10% of the density values for each 
bin, which are then scaled as per the profile into a mag- 
nitude density. Note that the slopes of the fits are nearly 
identical to those from analysis of the integrated-light pro- 
file, although the results are very slightly shallower. 

There are striking differences between the EFF fits 
to the MAD data and the optical results from |Mackey fc| 
~ (|2003|). We also fitted EFF profiles to the data 



Gilmore 

from |Mackey fc Gilmore] with it truncated to the outer 
radius of the MAD data. However, without measurements 
closer to the core it is difficult to obtain a meaningful com- 
parison for /io and a. Their inner data points suggest a 
somewhat steeper turnover, with |Mackey fc Gilmore| report- 
ing 7 = 2.43 ±0.09 for their inner EFF component, but the 
integrated-light MAD measurements at r < 1" rule out such 
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Figure 17. Azimuthal density profiles for a maximum radius of 
28'.'0 from the core of R136 and a completeness cut-off of 50%, 
demonstrating that the cluster is not relaxed over a ~7pc scale. 
Due north from the cluster centre corresponds to <f> = 0°. 
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Figure 18. As Figure |17| but with a completeness cut-off of 40%. 



a slope in the near-IR. Most relevantly, fits to the slope of 
the optical data yield 7 =1.80 ±0.10 in both the V and I 
bands, in reasonab le ag reement with the MAD data, as il- 



lustrated in Figure 21 ' Note the presence of the bump (or 



dent) noted by both |Meylan| ( |1993| ) and |Mackey fc Gilmore] 



(20031, at a radius of approximately 10 . For comparison, 
McLaughlin fc van der Marel| (20051 found 7 = 2.05 from 
fits to the full range of the HST dataQ 

A value of 7 = 1.65 ± 0.08 was also found for the slope of 
the number density profile for stars brighter than Mr = —5 
by Moffat et al. | ( |1994[ ) . The extension of the power-law out 
to ~ 10 pc makes R136 a good example of clusters with an 
extended halo - not all young clusters have such a halo, 



as illustrated by Moffat et al. ( 1994 1 through a comparison 



6 These are the original data from |Mackey fe Gilmore||McLaugh-1 
|lin &: van der Ma rcl (2005) later identified a magnitude offset of 
—0.715 (i.e. a brighter profile) compared to Mackey & Gilmore, 
w ho had masked out bright star s. 

7 |McLaughlin &: van der Marel| | |2005| find 73D =3.05, which is 
the slope of the 3D density profile, where 73D = 7 + 1 because 
of projection). 
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Figure 19. EFF fits to the combined MAD profiles for Field 1 
(with the iC s -band profile offset to fainter magnitudes for clarity). 
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Figure 20. EFF fits to the combined MAD profiles for Field 3 
(with the ii's-band profile offset to fainter magnitudes). 
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Figure 21. Profile fits to the HST optical data from Mackcy & 
Gilmorc (2003), with the /-band profile offset to fainter magni- 
tudes. 



with NGC 3603, which lacks bright stars beyond ~ 1 pc. 
|Mafz Apell aniz ( 200l| shows several extra-galactic examples 
of clusters with and without an extended halo. The origin 
of such halos in some clusters compared to those without 
remains unclear at present. 



6.4 Discussion 

6.4-1 Structural parameters of R136 

In summary, we find slopes of 7 < 2 from EFF fits to the 
near-IR luminosity profile of R136. If extrapolated to infinity 
such a profile would have infinite mass, so this can not extend 
to very large radii as the bulk of the mass would then be 
far from the core. It remains to be seen what the origin of 
these very shallow profiles is, although |Larsen| ( |2004| ) found 
that shallower profiles were more common in the younger 
members of their extragalactic cluster sample, and that 7 
appears to increase with age. The development of a core 
and the reduction of the central density with age is likely 
the result of mass loss by stellar evolution and heating by 
binaries in the centre ( jTakahashi fc Portegies Zwart|2000| ). 

Due to the very young age of R136 its density profile 
could be related to the formation process of the cluster, i.e. 
while it is still in its earliest stages of evolution we might see 
the imprint of the parent molecular cloud. Molecular clouds 
are approximately isothermal spheres in which, with 7 = 1 



(i.e. 730 = 2), the density scales as pocr 2 (e.g. McGlynn 
[1984] . For example, |Huff fc StahTer] ( [2006] ) found that the 
Orion Nebula Cluster (with an age of ~1 Myr) can be well 
approximated by a power-law profile with 7 = 1. Alterna- 
tively, the dissipationless collapse of a (gas free) stellar sys- 
tem also results in a power-law density profile, and possibly a 
core. The size of the core and the slope of the density profile 
depends mainly on the virial temperature in the initial con- 
figuration (McGlynn 198 4}. For a co ld collapse, i.e. zero ini- 
tial velocities of the stars, |McGlynn| found ~ 1.5 < 7 < ~ 2, 
with a core size of almost zero. The 'warmer' the initial con- 
ditions, the steeper the density profile and the larger the 
size of the core in the final configuration. So the observed 
density profile in R136 could also be the result of a cloud 
collapse, with star formation occurring in the early phases of 
the collapse, such that the dynamical interactions between 
the stars determine the properties of the cluster. 



Hunter et al. (1995) found no evidence for mass seg- 
regation in R136. 1 Andersen et al.| ( |2009[ ) also ruled it out 
for r>3pc (=12'.'4), finding no evidence for the flattening 
of the IMF below 2M Q reported bylSirianni et al.l (12000b. 



An early AO study o f one 12'.'8 x 12'.' 8 quadrant of R136 by 
|Brandl et al.| ( |1996| ) sampled a similar mass range to the 
MAD data, and found evidence of mass segregation above 
12 Mq. Whether the MAD data shows effects of the more 
massive stars being preferentially in the core will be ad- 
dressed elsewhere, but if evidence for mass segregation were 
found in the inner 3pc, then the core radius derived here 
would be somewhat smaller than the 'primordial' radius. 



6.4-2 Structure at larger radii 

Evidence for a second component in the near-IR profiles is 
less compelling than in the optical data, and we do not in- 
voke a second component in our EFF fits. There is tenta- 
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tive evidence in Field 1 for a 'dent' at comparable radius 
to the optical data (Figure [19] cf. Figure [2lj ), but the data 
for Field 3 are fit more cleanly by a single-component pro- 



file (Figure 201. Unfortunately the MAD data do not lend 



themselves to analysis of the integrated light out to larger 
radii. While the 'stitching together' of the two parts of the 
profile can influence the appearance of the combined pro- 
file, exclusion of one or two bins in the overlap region when 
combining them has very little influence on the EFF results. 
Indeed, now armed with knowledge of the minimal core, the 
EFF fits to the optical data are well matched at larger radii. 

The lack of a second component is interesting since the- 
oretical explanations for bumps (or halos) around young 
dense clusters (e.g. |Bastian fc Goo dwin 2006| should be in- 
dependent of the wavelength in which the profile is con- 
structed. The temptation is to attribute the past (two- 
component) results to variable extinction in the optical im- 
ages. From the optical image of R136 (Figure [1} there is a 
notable 'void' to the north-east (running from north-west to 
south-east) at r ~ 10", suggesting that differential extinction 
could be a prime suspect. We also note that more prosaic 
reasons (using circular annuli to determine the luminosity 
profiles in elliptical clusters) are now also being suggested 
to account for apparent breaks in density profiles in other 
clusters ( jPerina et aT]|2009 L 

In this context, the asymmetric structure seen in Fig- 
ure [17] must contribute strongly - almost certainly leading 
to the differences seen in the outer regions in Figures [19] 
and |20| Further support of this hypothesis is provided by 
inspection of the original Wide Field Planetary Camera 2 



(WFPC2) images used by Mackey & Gilmore R136 was lo- 
cated slightly off-centre in the Planetary Camera detector, 
with the Wide-Field Cameras primarily sampling the region 
east of R136, i.e. Field 1 observed with MAD. Interestingly, 
|Andersen et al.| ( |2009[ ) only fit their EFF profile to the inner 
2pc (=8. 25), noting in the caption to their Figure 12 that 
the presence of individual bright stars introduces 'the jitter 
in the surface brightness profile', hinting at the arguments 
that we advance here. 

Such deviations from symmetry are not surprising in 
such a massive, intricate, multi-population and, as yet, seem- 
ingly unrelaxed star-formation region. This point is rein- 
forced by the identification of (at least) five distinct pop- 
ulations in 30 Dor by |Walborn fc Blades] (1997} , three of 
which are sampled by the MAD observations: the central 
'Carina Phase' (rich in early O-type stars and including 
R136), an older 'Scorpius OBI Phase' of early-type super- 
giants throughout the central field, and a young (likely trig- 
gered) 'Orion Phase' to the north and west, mostly embed- 
ded in the gaseous filaments visible in Figure [l] and partially 
observed by Field 2. The stars associated with these popu- 
lations in the region immediately around R136 are neatly 
illustrated by Figures 3 and 4 from ? - their distribution is 
far from uniform, illustrative of the asymmetries shown in 
Figure [l7j 

The twin populations around R136 invite the question 
of their formation history. Is the Sco OBI Phase simply 
the remaining 'field' population, dispersed from a previous 
star-formation event? Are the younger, more massive Ca- 
rina Phase stars formed locally, or are they ejected members 
of R136? The answer to the latter question is most likely a 
combination of both. Without comprehensive stellar and gas 



dynamics of the different populations in 30 Dor we are lim- 
ited to speculation at the present time, but this is one of 
the principal motivations for the VLT-FLAMES Tarantula 
Survey ( |Evans et al.|2010 l, which has obtained multi-epoch 
optical spectroscopy for over 1,000 stars and will address 
these points. 

The formation history of these 'halo' stars goes to the 
heart of the theory of infant mortality, where deviations from 
EFF profiles in M82-F, NGC 1569-A and NGC 1705-1 were 



argued by Bastian & Goodwin ( 2006 1 to arise from rapid gas 
removal in clusters undergoing violent relaxation. Note that 



the HST profile from |Mackey fc Gilmore] ( |2003| > extends be- 
yond a radius of 60", thus also sampling the triggered gener- 
ation as well as the halo around R136, i.e. some of the excess 



light reported by Mackey & Gilmore (and, by inference, by 
Bastian fc Goodwin[ ) is likely attributable to ejected stars, 
but that triggered star formation also contributes. 



To echo the sentiment of Walborn & Blades (19971, 



the central 2' of 30 Dor observed with MAD corresponds 
to a physical size of ~30pc. If one projects this to a dis- 
tance of 10 Mpc, the same structures would all be contained 
within ~0'.'5. This reinforces the need to take potential clus- 
ter asymmetries and triggered star-formation into account 
when interpreting the radial luminosity profiles of distant 
unresolved clusters. 

The MAD data have provided us with a truly unique 
view of the central region of 30 Doradus. In the broader con- 
text, the delivery of impressive image quality and stability 
across wide fields (compared to classical AO observations) 
is very encouraging in advance of future multi-laser, MCAO 
systems on 8-m class telescopes, and the development of 
wide field, AO-corrected, imagers and spectrometers for the 
E-ELT. 
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